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NOMENCLATURE 

a  -  blade  section  lift  curve  slope,  per  radian 

A,C,  ■>-  aerodynamic  center 

C^Qcl  -"  lift  deficiency-  function,  including  blade  wake  effects 

C.G.  -  center  of  gravity 

c  -^  blade  section  chord 

E  ^  blade  bending  modulus  of  elasticity 

E.E,  ^  elastic  axis 

g,  -^  displacement  of  blade  k  bending  mode 

r  -  blade  section  second  moment  of  area,  flatwise  bending, 

If 
or  — 

h 

I,  -  blade  moment  of  inertia  about  flapping  axis 

I-  -  blade  section  moment  of  inertia  about  its  elastic  axis 

L  -  nondimensional  gereralized  mass  of  k  bending  mode 

I  -  blade  section  moment  of  inertia  about  its  center  of 

o 

gravity 
I  -  nondimensional  blade  product  of  inertia  about  flapping 

and  elastic  axes 

MCr)  -  blade  bending  moment  at  r 

M.^  -  blade  section  aerodynamic  moment  about  elastic  axis 

Mg,MQ,Mo       -  nondimensional  aerodynamic  feathering  moment  derivatives 

3*  g'  g 

mg,mg,mg       ~  nondimensional  aerodynamic  derivatives  associated  with 
in2»%2'™02        blade 


'"e2''"62'™g 


m 
g 
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m  -  blade  mass  per  unit  length 

R  "Ts  blade  radius 

r  '^  radial  station  along  blade 

T  -  blade  thrust 

U  -^  velocitj^  perpendicular  to  blade  span  in  vertical 

plane  at  blade  element 
U„  -^  velocity  perpendicular  to  blade  span  in  horizontal 

plane  at  blade  element 

X  ^  ^/r 

X.  •  distance  between  A.C,  and  E.A.,  positive  forward 

Xj.  -  distance  between  C,G.  and  E.A,,  positive  forward 

z  ~  vertical  displacement  of  E.A.  of  blade  element  from 

horizontal 
B  r>  displacement  of  blade,  rigid  flapping  mode 

Y  •-  Locke  number  — 

^b 

\Cx)  -  mode  shape  of  blade  k  bending  mode 

9  -  blade  pitch 

V  -  nondimensional  root  of  bending-torsion  characteristic 

equation 
V.  -  nondimensional  rotating,  natural  frequency  of  k  bending 

mode 
p  -  air  density 

Wi^  -  rotating  natural  frequency  of  k  bending  mode 

Wq  «  nonrotating  natural  frequency  of  blade  torsional  motion 
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SUMMARY  AND  INTRODUCTON 

This  report  presents  the  results  of  a  theoretical  approach  to  the 
analysis  of  the  flutter  of  extremely  flexible,  helicopter  rotor  blades. 
The  research  was  initiated  by  the  fact  that  during  experiments  conducted 
at  the  United  States  Naval  Academy,  1969-1970,  it  was  found  that  there 
were  no  theoretical  or  experimental  valves  available  for  comparison.  Some 
research  had  been  done  in  the  derivation  of  equations  of  motion  for  a 
flexible  rotor  blade,  but  actual  numerical  calculations  had  not  been  made. 
With  this  in  mind,  the  equations  of  motion,  incorporating  three  degrees  of 
freedom,  were  derived  for  a  rotating,  flexible  rotor  blade.  The  three  degrees 
of  freedom  used  were  pitching,  flapping  and  the  first  elastic  mode  of  flap- 
bending.  The  equations  were  then  expanded  into  the  characteristic  equation. 
The  objective  of  the  report  is  to  analyze  the  effect  of  the  modulus  of  elasticity 
upon  the  flutter  boundaries  of  a  rotor  blade.  Also,  at  each  modulus  of  elas- 
ticity, it  is  desired  to  study  the  effects  of  the  movement  of  the  aerodynamic 
center  and  the  center  of  gravity  upon  the  flutter  boundary. 

The  area  of  work  is  done  entirely  in  the  field  of  classical  flutter. 
It  is  assumed  that  the  inflow  of  the  blade  is  not  stalled.  Hence,  there  is 
no  wake  or  stall  flutter  to  consider. 

The  theoretical  results  consist  primarily  of  stability  plots  made  for 
each  different  modulus  of  elasticity.  These  plots  are  arrived  at  through  the 
use  of  Routh's  criteria.  All  of  the  results  obtained  are  in  very  close  agree- 
ment with  personal  observations  made  by  the  author  o^  experiments  conducted  at 
the  United  States  Naval  Academy  CReference  C23)).  They  also  agree  with  obser-- 
vations  made  by  the  Martin  Company  (Reference  C12)),  However,  the  results  do 
vary  somewhat  from  what  previous  research  has  shown  to  be  true  for  fairly  rigid 


blades,  thus  implying  that  more  work  should  be  done  in  the  area  of  extremely 
flexible  rotor  blades. 


DERIVATION  OF  EQUATIONS  OF  MOTION 
Since,  in  this  report,  it  was  desired  to  study  the  effects  of 
modulus  of  elasticity  on  a  rotating  helicopter  blade,  it  was  necessary 
to  derive  equations  of  motion  which  would  incorporate  this  as  a 
variable.  This  necessitated  having  the  equations  being  of  more  than 
two  degrees  of  freedom  in  order  to  incorporate  the  bending  effects 
of  the  blade.  Three  degrees  of  freedom  were  considered;   the  first 
torsional  mode,  the  zeroth  flapping  mode,  and  the  first  elastic  mode 
of  flap-bending.  The  two  flapping,  modes  are  illustrated  in  Figure  23. 
Ideally,  of  course,  an  infinite  number  of  degrees  of  freedom  should 
be  used  for  the  analysis,  but  the  mathematics  become  impossible  to 
handle  as  can  be  seen  in  Reference  Cl)>  which  only  used  ten  degrees 
of  freedom.  Also,  it  was  decided  to  use  only  three  degrees  of  freedom 
because  it  was  felt  that  these  were  the  ones  which  Imve  the  primary 
effect  on  blade  motion. 

The  method  of  the  derivation  of  the  equations  of  motion  vary 
from  reference  to  reference.   For  example,  References  (2),    C21),  and 
(11)  use  LaGranges  method  for  the  derivation,  whereas  References  CIS} 
and  (5)   used  the  virtual  work  method.  However,  the  method  used  in 
this  report  follows  that  of  References  (18)  and  (13)  and  is  the  sum- 
mation of  forces  or  moments. 

Prior  to  the  derivation  a  number  of  assumptions  concerning  the 
rotor  blade  were  made.  They  are: 

U)   The  rotor  blade  is  untwisted* 

C2)  Only  the  first  torsional  mode  is  considered, 

C3)  The  rotor  blade  is  untapered. 


(4)  The  rotor  blade  is  flexible,  or  in  other  words,  the  first 
elastic  mode  is  considered. 

(5)  The  rotor  blade  is  of  a  uniform  mass  distribution. 

(6)  The  rotor  blade  is  hinged  at  the  root,  or  is  articulated, 
with  no  hinge  offset. 

Further  assumptions,  as  needed,  will  be  made  as  the  derivation 
progresses. 

First,  consider  Figure  1,  which  shows  the  bending  out  of  the 
plane  of  rotation  of  a  flexible  rotor  blade.  The  positive  direction 
for  "z"  will  be  taken  as  a  bending  of  the  blade  upward.   It  is  ev- 
ident that  the  bending  moment  at  r  due  to  the  forces  at  s  will  be 

To  arrive  at  the  differential  equation  for  beam  bending,  equation 
Cl)  must  be  differentiated  twice  with  respect  to  r,  which  yields 


(2) 


d  K^       d  ^  d  1^      Jt- 


-  d^M  ^ri  m(h) 


d  r 
and  from  beam  theory 


Therefore,  the  differential  equation  for  a  blade  bending  out  of 
the  plane  of  rotation  in  a  centrifugal  force  field  is,  after  rearrange- 
ment 


d^' 


iJ_  1  m  S  fl^d  S 


(3) 


Th.e  form  that  equation  C3}  is  in  does  not  lend  itself  to  a  flutter 
analysis.  Hence,  with  this  in  mind,  a  series  solution,  written  in 
terms  of  normal  bending  modes  was  assumed: 


(4) 


For  the  moment,  assume  that  the  aerodynamic  damping  is  zero  —  i.e. 
dT(r)/dr  =  0.  Then,  substituting  equation  (4\   into  equation  (3)   the 
bending  equation  becomes 


oO 


<=jL 


c/  *-  V    d  ^ 


4-  r>r)n^d_K^J  g^  f  VI  .^/g;  I 


C5) 


-  o 


assLune  that  the  motion  is  simple  harmonic  at  the  k   rotating,  undamped 
natual  bending  frequency.  Therefore: 


gi^^gK 


iO,Xl.1- 


■2   a      /■^^k^^       9   2 


hence   [   J  gk  -  ^k  ^^  '^1'<  S*^ 


Now,  putting  the  aerodynamic  damping  back  into  the  equation,  and 
with,  the  above  assumption,  the  bending  equation  of  motion  becomes 


'\^^'^K^\%   ^   l^^^-^X^   = 


It  is  now  necessary  to  consider  the  effect  of  torsional  motion 
upon  the  bending  equation.  There  is  an  inertial  term  and  a  centrifugal 
force  term  due  to  the  torsional  motion  which  causes  the  coupled  blade 
bending  -  torsion  equation  of  motion  to  be 


.oO  . oO 


Z^.i-^nK^t  +  Z,..-L^^*^''-^^^k3K  ->^h'^ 


(7) 


multiplying  by  x\.    and  integrating  from  0  to  R  we  get 

a 


where 


The  value  of  k.  is  dependent  upon  the  number  of  modes  of  bending 
chosen.  In  this  report  k  =  2  since  the  rigid  flapping  and  the  first 
elastic  mode  are  the  two  modes  of  bending  being  considered. 

It  will  be  assumed  that  the  bending  motion  can  be  described  by 

^=•-13^  19  ''' 

where         r  =  rigid  flapping  mode 

T\   =  first  elastic  mode 
By  substituting  (9)    into  (8)   and  nondimensionalizing  the  resulting 
two  equations  by  dividing  by  I^^^  we  get 

I^e  r  J-  \   t-Uji-  (10) 


i'-  ^^       "n"      ■     -  LA^j„  d^ 


h"^ 


where  X^^tAx 


X.- 


lb 


and  also,  the  nondimensional  rotating  natural  frequency  of  the  rigid 
flapping  mode,  v-^,    is  equal  to  1.0, 

Equations  (10)  and  (H)  are  the  two  bending  equations  of  motion 
for  a  rotating,  flexible  rotor  blade. 

Next,  the  blade  torsional  equation  of  motion  must  be  derived. 
Figure  2  shows  the  torsional  geometry  used.   It  can  be  seen,  from  the 
figure,  that  the  positive  direction  for  "e»«  has  been  chosen  as  a  pitch 


dovm  of  thje  blade. 

TaleLng  moments  about  the  elastic  axis 


-Jo     cir 

2 
substituting  for  z  and  nondimensionalizing  by  I,fi  we  get 


013) 


Equation  Q^),  then,  is  the  torsional  equation  of  motion  for  a 
rotating,  flexible  rotor  blade. 

It  is  now  necessary  that  the  aerodynamic  terms  of  equations  CIO), 
Cll),  and  C13)  be  expanded.  From  Reference  C3),  page  272,  for  a  unit 
span  of  a  thin  airfoil  oscillating  in  an  incompressible  flow  we  have 


-  1,  ja.o^T-  C'W[Up^(>!^o  t(.^c-A;,)e] 
2 


C14) 


d^-     8 
-3^'  e-(^A-.E^c)V J  +  f  J«-c  C^ A-^  CM-  CIS) 

where 

The  term,  C'O^),  is  a  modified  Theodorsen  function  Csee  Reference 
(16)).  It  has  been  modified  to  take  into  account  the  wake  effect  of 
the  rotating  rotor  blade.   In  this  report,  it  is  assumed  that  C'Ck) 
has  a  value  of  one.  Since  this  report  is  concerned  only  with  the 
hovering  case,  this  assumption  is  a  very  good  approximation. 

Normally,  in  this  type  of  analysis,  the  mass  terms  which  are  pro- 
portional  to  z  and  9  are  neglected,  as  they  shall  be  in  this  report. 
If  it  is  desired,  the  effect  of  these  terms  can  be  included  by  the  ap- 
propriate modification  of  the  values  of  the  blade  inertial  constants, 
but  the  change  is  usually  of  a  negligible  magnitude. 

From  References   (18)  and  (21),  the  aerodynamic  terms  can  be 
expressed  as 


where    /^  =  8    CM   [it  JyU  ■;(»/  f  J 

Tit)  [C 1  V'-'  +  f/a  S  -N  V  -^  2  cost2Vj7 

=  f  cYtF  [i>  COS  f  +  ^^  "=  s  'N  C2«] 


n^ 


^.  =Z  ^W  [CCtf^dx  4^siN  V  [;ct)^cix 


where 

^,,=  f  c"W  [  i'ai^^dx  f^^si^^r  j^^^JArJ 
f    J!  c'toy*  S'W'^X    'C^d'^' 


C17) 
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m„=  f  C'(K)[^c.sViM  ^©J'*'  ^T/<'sivU'«{tfx(^)J^] 


and  finally. 


where 

Me  =  8- [r  -  f  Cte   ^  J  [^  -  ^]  [t  V  Si*"^J 

M6=  -6-1^  C^[l+  3yts/A/Vt  S^^S/a;>] 


C18) 


u 


Although,  these  aerodynamic  expressions  look  very'  complicated, 
by  assuming  that  the  rotor  is  operating  in  the  hovering  state,  y 
equals  zero.  This  fact  greatly  simplifies  the  above  expressions. 
Also,  if  c«R,  then  M^  and  M*  are  normally  set  equal  to  zero. 

By  assuming  that  the  advance  ratio  is  zero,  the  following 
derivatives  become  zero: 
M3  =  0 

Mg2  =  0 

M32  =  ° 

M  =  0 
g 

M3  =  0 

M  =  0 
g 

Hence,  the  final  form  of  the  equations  of  motion  become:  the  flapping 

equation 


CI  9) 


the  flap-bending  equation 


■^^\   ^  ^^4-  +  ^^^^'''3  ~  ^xi^^  ^  Ki-i'xO® 


C20) 


0- 


12 


and  the  torsion  equation 


-Ix|_^M^|.-J,^^i|.,  ^M,i. 


+  Ktr+itt)^]e-L,4_fM^^-r^9--o     ^''- 


O-^  ^  A 


This  concludes  the  derivation  of  the  equations  of  motion  for  a 
flexible,  rotating  rotor  blade.  Appendix  IV  gives  an  explanation  of 
the  various  terms  in  the  above  equations. 
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ANALYSIS  OF  EQUATIONS  OF  MOTION 

There  are  a  couple  of  ways  in  which  the  equations  of  motion  can  be 
analyzed  in  order  to  obtain  the  flutter  boundaries.   By  one  method,  the 
equations  of  motion  have  to  be  solved  for  the  critical  flutter  frequency 
and  the  RPM  at  which  the  flutter  instability  begins.  This  is  a  very 
good  approach  to  the  problem  as  long  as  there  are  no  more  than  two  degrees 
of  freedom.  As  soon  as  the  system  is  increased  by  even  only  one  degree 
of  freedom,  the  mathematics  involved  in  the  solution  become  unmanageable, 
without  the  aid  of  a  lot  of  computer  time.  Due  to  the  time  restriction 
placed  on  this  report  this  method  was  ruled  out  as  just  being  too  time 
consuming.   If  there  is  interest  in  this  approach.  References  (8)   and 
C19)  would  be  of  great  advantage, 

A  second,  much  simpler  method  of  obtaining  the  flutter  boundaries  is 
to  use  Routh's  Criteria.  This,  as  in  the  first  method,  also  requires  that 
the  characteristic  equation  be  obtained.  However,  a  solution  of  the  equa- 
tion is  not  needed.  According  to  Routh,  for  a  system  to  be  dynamically 
stable  it  is  necessary  that  all  the  coefficients  of  the  characteristic 
equation  be  of  the  same  sign.  Also,  at  the  same  time,  it  is  necessary 
that  Routh's  discriminant  also  be  positive.  This  discriminant  is  obtained 
from  an  array.  The  development  of  this  array  and  the  discriminants  is 
shown  in  Appendix  II. 

Now  that  the  method  of  analysis  is  known,  the  equations  of  motion 
must  be  put  in  the  correct  form.  In  other  words,  the  characteristic  equa- 
tion of  motion  has  to  be  obtained.   Before  this  is  done,  though,  it  should 
be  pointed  out  that  there  are  no  time-dependent  terms  in  the  equations  of 
motion.  This  is  because  the  hovering  condition  has  been  assumed  and  hence 
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the  advance  ratio  is  zero  which,  eliminates  the  time-dependent  terms.  This 
means  that  all  of  coefficients  of  the  characteristic  equation  will  be 
constant,  and  that  a  conventional  technique  for  an  explicit  solution  of 
the  equation  could  be  used  if  so  desired.  Even  though  solution  was  not 
attempted,  having  constant  coefficients  simplified  the  analysis. 

To  further  simplify  the  analysis,  it  was  assumed  that  the  motion  of 
the  blade  was  a  simple  harmonic  motion.  Therefore,  it  was  known  that  the 
harmonic  perturbations  of  the  variables  B,  9,  and  g  would  be  of  the 
following  form: 

S  '  ^o  e  C22) 

3=  ^°^  C23) 

©  ^  ©0^  C24) 

Substituting  (22),  (23),  and  (24)  into  (19),  (20),  and  C21),  and  after 
dividing  all  the  equations  by  e   ,  the  equations  of  motion  become: 

the  flapping  equation 

the  flap  -  bending  equation 


M^"^    -f^n^-^V  ^I^V^^,  'I^^e.-^^^ 


(26) 
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and  the  torsion  equation 

It  can  be  seen  that  now  all  three  equations  are  in  terms  of  v, 
which  is  a  non-dimensional  root. 

V  =  w 

If  the  equations  were  solved  for  v,  the  answer  would  be  the  critical, 
reduced  frequency  for  flutter.   Since  the  equation  was  not  to  be  solved  an- 
other parameter  had  to  be  chosen  for  the  stability  plots.  This  was  the 
nonrotating,  reduced  torsional  frequency,   /fi. 

The  equations  can  now  be  put  into  a  determinant  form,  and  the  de- 
terminant expanded  to  arrive  at  the  characteristic  equation.  This  also 
shows  why  more  than  a  few  degrees  of  freedom  are  not  very  often  analyzed. 
With  three  degrees  of  freedom  the  determinant  is  a  3  X  3,  and  the  expansion 
is  fairly  complicated.  As  the  degrees  of  freedom  are  increased  the  square 
determinant  increases  by  one  and  becomes  more  and  more  difficult  to  expand. 
The  full  expansion  of  the  determinant  used  in  this  report  may  be  found  in 
Appendix  I.  A  quick  look  at  Appendix  I  will  show  how  unwieldy  the  coef- 
ficients of  the  characteristic  equation  become,  with  some  coefficients 
having  as  many  as  thirty  (30)  terms.   Increasing  the  system  by  only  one  de- 
gree of  freedom  nearly  doubles  the  number  of  terms  of  the  coefficients,  and 
the  complexity  of  the  problem. 

After  expansion  of  the  determinant,  the  characteristic  equation  is 
found  to  be  a  sixth-order  equation  of  the  following  form: 
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/\\?S^v7^^  C^^'^-h  D^'^t-  E\?^  i-  ¥^   -^G-C  (28) 

Once  the  characteristic  equation  was  derived,  it  was  possible  to 
calculate  the  coefficients  for  the  varying  parameters  and  to  establish  the 
flutter  boundaries.  The  main  parameter  of  interest  was  the  modulus  of 
elasticity  of  the  rotor  blade.  Since  the  equations  of  motion  were  derived 
in  terms  of  normalized  mode  shapes,  the  only  place  the  modulus  of  elas- 
ticity enters  the  equations  is  through  the  mode  shapes.  Hence,  it  was 
necessary  to  use  several  different  mode  shapes  in  order  to  study  the  ef- 
fect of  a  modulus  of  elasticity  tending  towards  zero.  The  specific  mode 
shapes  used  and  why,  is  explained  elsewhere  in  this  report. 

Besides  studying  the  effect  of  the  modulus  of  elasticity,  it  was  also 
desired  to  study  the  effects  of  other  parameters  at  each  respective  mode 
shape.  At  each  mode  shape,  two  other  parameters  were  varied.  These  were 
the  center  of  gravity  location  and  the  aerodynamic  center  location.  The 
location  of  both  of  these  axes  was  with  respect  to  the  elastic  center  axis. 
l\\e   center  of  gravity  and  the  aerodynamic  center  were  varied  first  inde- 
pendently and  then  they  were  varied  together.  All  of  the  results  were  then 
plotted.  In  all  the  plots  the  ordinate  is   /f2,  which  is  the  reduced,  non- 
rotating  natural  frequency  of  blade  torsional  motion.   By  using  this  as  a 
parameter  it  is  possible  to  tell,  for  a  particular  torsional  frequency,  at 
what  RPM  the  blade  system  becomes  unstable.  Also,  for  a  particular  RPM  it 
is  possible  to  calculate  the  corresponding  critical  nonrotating  torsional 
frequency.  The  abscissa  of  the  stability  plots  was  either  the  center  of 
gravity  location  or  the  aerodynamic  center  location. 
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Due  to  the  complexit)^  of  the  coefficients  of  the  characteristic 
equation,  all  the  numerical  calculations  were  done  on  a  computer. 
The  computer  used  was  IBM  36(1  Model  91.  The  computer  was  located  at 
the  Princeton  University  Computer  Center  at  Princeton  University.  Since 
the  program  used  was  just  a  simple  mathematical  routine,  it  is  not  in- 
cluded in  this  report.  All  of  the  results  were  transferred  to  graphs, 
which  are  included,  however. 
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MODE  SHAPES 

As  lias  already  been  seen,  the  equations  of  motion  of  the  rotating 
rotor  blade  were  derived  in  terms  of  normalized  mode  shapes.   It  should 
also  be  noticed  tfiat  the  modulus  of  elasticity  does  not  appear  explicitly 
in  the  equation,  but  rather  is  incorporated  into  the  mode  shapes.  This 
made  it  necessary  to  choose  various  mode  shapes  which  would  be  representa- 
tive of  the  varying  modulus  of  elasticity.   It  was  possible  to  obtain  ex- 
act mode  shapes  for  two  limiting  cases  Qsee  References  QS}  and  (l^)) . 

The  two  limiting  cases  are 

(XX   A  rotor  blade  with  the  modulus  of  elasticity  equal  to  zero. 

C21  A  rotor  blade  that  is  considered  to  be  a  typical  rigid  rotor  blade. 

It  was  not  possible  to  find  exact  mode  shapes  for  a  modulus  of  elas- 
ticity between  case  (X)   and  case  C2) .  However,  by  graphing  the  above 
two  mode  shapes  and  then  drawing  curves  which  fell  between  these  two  cases, 
it  was  possible  to  obtain  other  mode  shapes.  See  Figure  3.   Only 
two  intermediary  mode  were  chosen,  making  a  total  of  four  (A) . 

The  reason  that  an  infinitely  stiff  blade  was  not  chosen  as  the  lim- 
iting case  was  that  it  is  nearly  impossible  to  achieve  such  a  situation 
with  a  long,  thin  rotor  blade.  Therefore,  the  mode  shape  of  a  typical 
rotor  blade  was  used.  Mathematically  the  mode  shape  can  be  expressed  as: 

^=  3X-  4-X^  C29) 

This  mode  shape  was  then  substituted  into  the  nondimensional  derivatives 
giving  the  following: 


?r,^-6 


M,-li[i-^^]fe-^J 
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These  derivatives  were  then  used  to  calculate  the  coefficients  of  the 
characteristic  equation. 

Now  look  at  case  (l\\     the  rotor  blade  with  the  modulus  of  elas- 
ticity equal  to  zero.  The  mode  shape  for  this  case  was  evaluated  in  a 
Ph.  D.  thesis  written  at  Rensselaer  Polytechnic  Institute.   CReference 
C14)).   It  was  found  that  the  mode  shape  was  of  the  following  form: 

|=^[3x-f)f']  C30) 

With  this  as  the  mode  shape  the  non-dimensional  derivatives  which  are 

affected  by  the  mode  shape  now  have  the  following  values: 
jY)   -  —  ~ — -  /m  -      '  ' —- 

% 
The  rest  of  the  derivatives  do  not  change  in  value. 

After  much  drawing  of  curves  and  curve  fitting,  it  was  possible  to  - 
arrive  at  two  other  mode  shapes  which  fell  between  the  two  limiting  cases. 
The  first  of  these  intermediary  shape  was  of  the  form: 

C3n 


t=  3X-  fTA^ 
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with,  the  corresponding  derivatives  having  the  values: 

The  second  intermediar)^  shape  was  of  the  form; 

\,f[^^-^K^l  C32) 

and  its  corresponding  derivatives  had  the  values 

One  interesting  fact  discovered  while  arriving  at  the  mode  shapes 
was  that  V-,  the  nondimensional  rotating  natural  frequency  of  the  second 
bending  mode,  did  not  change  as  the  modulus  of  elasticity  changed.   Csee 
References  0-^\   and  C20)}.  Hence,  it  was  possible  to  use  the  same  fre- 
quency throughout  all  the  calculations.   In  this  case  the  frequency  turned 
out  to  be  2.5. 

The  reason  there  is  no  change  in  the  frequency  is  that  it  is  more 
strongly  dependent  upon  mass  distribution  than  it  is  modulus  of  elasticity. 
And  since  all  the  cases  in  this  report  have  the  same  mass  distribution, 
i.e.  uniform,  the  frequency  of  the  second  bending  mode,  or  in  other  words, 
the  first  elastic  mode,  remains  constant. 

It  should  also  be  mentioned  that  the  normalized  mode  shapes  are  also 
strongly  dependent  upon  mass  distribution.  Although  this  fact  was  not 
considered  in  this  report,  it  would  be  worthwhile  to  do  so. 
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RESULTS 
Overall  Results 

As  we  previously^  mentioned,  the  coefficients  of  the  characteristic 
equation  and  the  values  of  Routh's  discriminants  were  all  calculated;  and 
then,  using  Routh's  criteria,  the  flutter  boundaries  were  plotted.   For 
all  the  calculations,  one  rotor  blade  type  was  used.  The  rotor  blade  used 
was  a  model  rotor  blade  developed  by  Cornell  Aeronautical  Laboratories  for 
some  of  their  experiments.   Csee  References  (6}   and  (S)) .     The  represent- 
ative data  for  this  model  blade  is  given  in  Appendix  III. 

For  each  mode  shape,  four  different  situations  were  analyzed  in  order 
to  provide  a  complete  flutter  analysis  of  the  rotor  blade  for  the  para- 
meters involved.  The  cases  investigated  were: 

CI)  The  aerodynamic  center  and  the  center  of  gravity  were  coincident 
with  the  elastic  axis. 

(2)  The  aerodynamic  center  was  coincident  with  the  elastic  axis,  and 
the  center  of  gravity  was  varied  in  location  forward  and  aft  of  the  elastic 
axis. 

C3)  The  center  of  gravity  was  coincident  with  the  elastic  axis,  and 
the  aerodynamic  center  was  varied  in  location  forward  and  aft  of , the  elastic 
axis. 

(4)  And  finally,  both  the  center  of  gravity  and  the  aerodynamic  center 
were  varied  in  location  forward  and  aft  of  the  elastic  axis. 
Each  case  shall  be  discussed  separately,  and  then,  an  overall  analysis 
shall  be  given. 

First,  let's  consider  case  Cl)»  where  all  three  axes  are  coincident. 
It  was  not  possible  to  plot  a  graph  of  the  dynamic  stability  of  this  system 


22 


in  a  manner  which  could  easily  he  understood.  Therefore,  the  results 
were  put  in  tahle  form  and  are  shown  in  Table  1.   In  the  table,  the  mode 
shapes  are  numbered  1  ttirough  4,  with,  mode  shape  number  1  being  the  most 
rigid  case  and  mode  shape  number  4  being  the  most  flexible  case.  The 
first  tiling  that  will  be  noticed  is  that  as  the  modulus  of  elasticity  is 
decreased  Cthe  blade  becomes  more  flexible)  the  RPM  at  which  the  rotor 
blade  becomes  unstable  is  decreased  Can  increase  in  /9.     is  a  decrease 
in  RPM) . 

From  Reference  (22)  it  was  found  that  a  typical  operating  value  for 
/9.   is  8.0.  Hence,  it  immediately  becomes  evident,  that  with  an  ex- 
tremely flexible  rotor  blade  the  blade  would  always  be  just  on  the  verge 
of  instability,  requiring  only  a  slight  abnormal  operating  condition  to 
cause  an  unstable  flutter.  To  make  matters  even  worse,  there  were  some 
observations  made  by  the  Martin  Company,  of  extremely  flexible  blades 
(Reference  (12))  which  showed  that  as  the  modulus  of  elasticity  of  a 
rotor  blade  is  decreased,  the  nonrotating  natural  frequency  of  the  blade 
torsional  motion,  0)^,  also  decreases.  This  means  that  /^   will  also  be 
reduced  and  hence,  the  stability  margin  is  decreased.  All  of  this  points 
to  the  problems  designers  have  with  extremely  flexible  rotor  blades. 

In  the  second  case,  the  aerodynamic  center  and  the  elastic  axis  were 
coincident.  The  results  are  shown  in  Figure  4.  Once  again  it  was  ob- 
served that  as  the  modulus  of  elasticity  was  decreased,  the  region  of  in- 
stability was  increased.   It  can  be  seen  that  there  was  a  gradual  increase 
in  instability  for  all  mode  shapes  as  the  center  of  gravity  was  displaced 
further  and  further  from  the  elastic  cixis.  However,  when  the  center  of 
gravity  was  located  more  than  20%  chord  from  the  elastic  axis  the  system 
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became  unstable  for  all  values  of  fj—  analyzed.  Other  References 
do  not  show  this  marked  instability  due  to  the  fact  that  in  the  majority 
of  the  analyses  the  center  of  gravity  is  not  moved  as  far  as  20%  chord 
away  from  the  elastic  axis.  The  reason  for  this  sharp  instability  is 
discussed  elsewhere  in  this  report. 

In.  case  (5),   the  center  of  gravity  was  coincident  with  the  elastic 
axis.  This  says  that  the  rotor  blade  was  mass  balanced.  Again  it  was 
observed  that  a  reduction  in  the  stiffness  of  a  rotor  blade  increases  the 
region  of  instability.   It  was  also  found  that  as  the  aerodynamic  center 
was  moved  aft  of  the  elastic  axis,  the  flutter  boundary  began  to  fall  off, 
whereas  if  it  was  moved  forward  of  the  elastic  axis  a  plateau  for  the  flutter 
boundary  was  reached.  See  Figure  5.  This  would  suggest  that  it  is  pre- 
ferable to  have  the  aerodynamic  center  aft  of  the  elastic  axis.  That  is 
the  same  conclusion  reached  by  Goldman  in  Reference  (12).  Contrary  to 
case  C2),  and  as  will  presently  be  seen,  case  (4),  case  (3)  does  not  have 
the  sharp  upward  break  in  the  flutter  boundary.  As  might  be  suspected, 
this  was  because  the  rotor  blade  was  mass  balanced. 

Figures  6-19  show  the  final  results  of  case  (4)j,  in  which  both  the 
center  of  gravity  and  the  aerodynamic  center  were  varied.  As  was  expected, 
it  was  found  that  in  this  case,  also,  that  reducing  the  stiffness  of  a 
blade  increased  the  instability.   It  was  further  found,  again,  that  it  was 
preferable  to  have  the  aerodynamic  center  located  aft  of  the  elastic  axis, 
(see  Figure  20) .  This  was  the  same  conclusion  reached  when  case  (3)  was 
analyzed.  Case  (4)  also  illustrated  the  same  strange  break  in  the  flutter 
boundary  that  was  found  in  case  (2) .   For  some  reason  there  appeared  to  be 
a  limiting  center  of  gravity  position,  beyond  which  the  system. was  unstable 
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for  all  values  of    /^. 

Since  an  optimum  positon  for  the  aerodynamic  center  had  been  found, 
it  was  deemed  desirable  to  also  find  an  optimum  position  for  the  center  of 
gravity.  With  this  in  mind.  Figure  21  was  drawn.   From  this  figure  it  can 
be  seen  that  the  region  of  instability  will  be  the  smallest  when  the  center 
of  gravity  is  coincident  with  the  elastic  axis.   By  looking  at  the  equa- 
tions of  motion  it  can  be  seen  that  in  this  situation  the  torsional  flutter 
is  uncoupled  from  the  flapping  and  flap-bending  flutter.   Csee  Reference 
(19)).  As  soon  as  the  rotor  blade  is  no  longer  mass  balanced  the  pitching 
and  the  flapping  of  the  rotor  blade  become  coupled  and  the  instability  is 
increased.   It  was  also  observed  that  a  forward  center  of  gravity  move- 
ment was  preferred  over  a  rearward  movement.  This  fact  is  further  supported 
by  observations  made  by  Goldman  in  Reference  (12). 

In  no  other  references  is  there  any  mention  of  forward  center  of 
gravity  movement  made.  And  if  the  rotor  blade  is  fairly  rigid  there  is  no 
real  reason  to  consider  this  fact.  However,  if  a  rotor  blade  being  designed 
is  extremely  flexible,  the  center  of  gravity  location  could  be  the  source  of 
many  problems. 

From  personal  observations  made  by  the  author  of  some  experiments  con- 
ducted at  the  United  States  Naval  Academy  with  extremely  flexible  blades. 
Reference  (23) »  the  critical  importance  of  center  of  gravity  location  was  seen. 
The  rotor  blades  used  had  originally  been  constructed  with  the  idea  of  having 
the  center  of  gravity  and  the  elastic  axis  coincident.  However,  it  was  found 
that  actually  the  center  of  gravity  was  aft  of  the  elastic  axis.  At  the 
time  this  was  not  of  too  much  concern.  However,  upon  conducting  the  experi- 
ments it  was  not  possible  to  keep  the  blade  stable.  This  theoretical  analysis 
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real  reason  to  consider  this  fact.  However,  if  a  rotor  blade  being  designed 
is  extremely  flexible,  the  center  of  gravity  location  could  be  the  source  of 
many  problems. 

From  personal  observations  made  by  the  author  of  some  experiments  con- 
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Reference  (23) »  the  critical  importance  of  center  of  gravity  location  was  seen. 
The  rotor  blades  used  had  originally  been  constructed  with  the  idea  of  having 
the  center  of  gravity  and  the  elastic  axis  coincident.  However,  it  was  found 
that  actually  the  center  of  gravity  was  aft  of  the  elastic  axis.  At  the 
time  this  was  not  of  too  much  concern.  However,  upon  conducting  the  experi- 
ments it  was  not  possible  to  keep  the  blade  stable.  This  theoretical  analysis 
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now  explains  why  tliis  was  so.  For  flexible,  i.e.  a  very  small  modulus 
of  elasticity,  the  center  of  gravity  should  be  coincident  with  or  for- 
ward of  the  elastic  axis. 

All  of  the  above  results  can  be  reduced  down  to  a  few  certain  facts 
that  must  be  kept  in  mind  when  designing  an  extremely  flexible  rotor  blade. 
Number  one  is  that  it  is  desirable  to  Fiave  the  aerodynamic  center  aft  of 
the  elastic  axis.  Secondly,  it  was  found  that  the  center  of  gravity  lo- 
cation is  critical  and  should  be  conicident  with  or  forward  of  the  elastic 
axis.   It  was  also  found  that  as  the  modulus  of  elasticity  is  reduced,  the 
nonrotating  natural  frequency  of  blade  torsional  motion,  o)  ,  becomes  an 
important  factor.  The  reduction  of  w  that  occurs  with  a  reduction  of 
modulus  of  elasticity  tends  to  place  the  operating  mode  of  the  rotor  blade 
a  lot  closer  to  instability.  This  also  explains  why  it  is  desirable  to 
have  the  center  of  gravity  forward  of  the  elastic  axis.  The  extra  margin 
of  stability  that  is  gained  is  needed. 

One  other  observation  made  in  this  analysis  concerned  the  effect  of 
center  of  gravity  displacement  from  the  elastic  axis  upon  the  torsional 
frequency.   From  the  figures  it  can  easily  be  seen  that  as  the  center  of 
gravity  is  displaced  further  and  further  from  the  elastic  axis,  in  either 
direction,  the  torsional  frequency  required  for  flutter  is  increased. 
Flutter  Phenomena 

As  was  found  in  discussing  the  overall  results,  case  (2)  and  case  (4) 
of  the  analysis  disclosed  a  very  sharp  upward  break  in  the  flutter  boundary 
that  occured  at  a  center  of  gravity  location  of  more  than  20%  chord  dis- 
placement from  the  elastic  axis.  Since  this  break  in  the  flutter  curve 
occured  within  one  chord  length  it  was  considered  necessary  to  analyze  it 
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and  see  if  the  reason,  for  its  occurrence  could  he  discovered.  First,  it 
was  important  to  find  out  whether  it  was  the  coefficients  of  the  charact- 
eristic equation  or  one  of  Routh's  discriminants  which  was  displaying  the 
instability. 

After  going  back  over  all  of  the  data,  it  was  found  that  in  all 
instances  it  was  the  coefficients  of  the  characteristic  equation  which  were 
displaying  the  instability.  Also,  it  was  observed  that  it  was  only  the 
coefficients  of  the  two  higFiest  order  terms,  or  in  other  words  the  "A"  co- 
efficient and  the  "B"  coefficient  if  the  equation  is  written  in  the  form  of 
Equation  (28).  The  next  step  in  the  analysis  was  then  to  find  out  which 
term  of  the  coefficient  caused  the  instability,  or  speaking  mathematically, 
which  term  caused  the  coefficient  to  become  negative. 

The  two  coefficients  may  be  expressed  as  follows: 

A=ir^-x,\-x;r,  C33) 

Although  all  cases  were  analyzed,  the  case  that  will  be  discussed  in 
this  report  will  be  the  case  where  the  modulus  of  elasticity  equals  zero. 
The  results  are  the  same  for  the  rest  of  the  cases,  but  with  a  little  dif- 
ferent numbers  for  the  coefficients. 

Using  the  derivatives  for  the  case  of  zero  modulus  of  elasticity,  and 
substituting  in  the  corresponding  numbers,  the  coefficients  may  now  be 
written  in  the  following  form: 


A^  ,ooozSl  '  .Ol^GXl  -  ,oB4-4-X^ 


(35) 
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i3-  -,01^  ^r  ^  .OZ-iX^X^~rniX^  f  ,0001  -.Olb'X;, 

+  ,Oh^X^^   f  .O001>qG-  ,OA-'hX^   -/-  ,0000  87  C36) 

Let*s  consider  the  "A"  coefficient  first.  A  look  at  the  numerical 
equation  shows  that  it  will  be  due  to  the  second  and  third  terms  that  the 
equation  first  becomes  negative  and  hence  the  system  becomes  unstable. 
These  two  terms,  wliich  are  dependent  upon  center  of  gravity  position,  are 
r  and  I  ^.   Both  of  these  terms  are  the  nondimensional  blade  product  of 
inertia  about  the  flapping  and  elastic  hinges.  Going  back  to  the  original 
equations  of  motion,  it  was  found  that  it  was  the  inertial  moments  about 
the  pitching  and  flapping  axes  which  gave  rise  to  the  above  terms.  These 
terms  are  also  called  the  mass  unbalance  terms  since  if  the  rotor  blade  is 
mass  balanced  they  are  zero. 

Now  let's  look  at  the  "B"  coefficient.   It  can  be  seen  that  this  co- 
efficient is  dependent  upon  both  center  of  gravity  position  and  aero- 
dynamic center  position.  However,  it  was  found  that  the  only  time  that  the 
"B"  coefficient  was  important  was  when  the  aerodynamic  center  was  aft  of 
the  elastic  axis  and  the  center  of  gravity  forward.  This  meant  that  the 
secomd  term  was  the  most  important  term.   It  was  also  found  that  I  had  the 
primary  effect,  as  it  did  in  the  "A"  coefficient. 

As  can  be  seen,  in  all  cases  the  sharp  break  in  the  flutter  curve  was 
due  to  the  mass  unbalance  terms,  I  and  I  -.   From  looking  at  the  equations 
of  motion  it  can  be  seen  that  the  instability  occurs  due  to  a  coupling  of 
the  inertias  between  the  pitching  and  flapping  of  the  rotor  blade.  Also, 
it  should  be  mentioned  that  these  two  inertia  moments  arise  due  to  either 
a  pitching  acceleration  or  a  flapping  acceleration. 

That,  then,  explains  the  reason  for  the  sharp  upward  break  in  the 
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All  this  shows  that  a  divergent  instability  becomes  of  no  more 
importance  for  an  extremely  flexible  blade  than  it  was  for  a  fairly 
rigid  rotor  blade. 
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CONCLUSIONS 
Based  on  the  theoretical  results  previously  discussed,  the  following 
basic  conclusions  were  reached. 

1.  As  the  modulus  of  elasticity-  of  a  rotor  blade  is  decreased  the 
region  of  instability  increases.  Since,  at  the  same  time,  w  is  decreasing 
there  can  be  a  very  real  problem  with  the  dynamic  stability  of  an  extremely 
flexible  rotor  blade. 

2.  In  the  design  of  an  extremely  flexible  rotor  blade,  it  is  desir- 
able to  have  the  aerodynamic  center  located  aft  of  the  elastic  axis  in  order 
to  obtain  dynamic  stability. 

3.  Also,  when  designing  an  extremely  flexible  rotor  blade  it  is 
necessary  to  have  the  center  of  gravity  coincident  with  the  elastic  axis 
or  forward  of  it.  This  is  contrary  to  classical  flutter  analysis,  but  is 
a  necessary  condition  if  the  blade  is  to  be  extremely  flexible,  due  to  the 
decrease  in  U)  . 

0 

4.  Coupling  of  the  mass  unbalance  terms,  or  in  other  words,  the  iner- 
tial  effects  due  to  pitching  and  flapping  acceleration  causes  a  sharp  up- 
ward break  in  the  flutter  boundary  if  the  center  of  gravity  is  displaced 
more  than  20%  chord  away  from  the  elastic  axis.  This  occurs  whether  the 
center  of  gravity  is  forward  or  aft  of  the  elastic  axis. 

5.  A  change  in  the  modulus  of  elasticity  of  a  rotor  blade  has  no 
effect  on  the  divergence  boundary  of  the  blade. 
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RECOMMENDATIONS 
Due  to  the  time  restriction  on  tlvis  report  there  were  many  areas  of 
the  flutter  problem  that  the  author  was  not  able  to  look  into.  With  this 
in  mind,  the  following  are  several  recommendations  for  future  study  in 
this  field. 

1.  Since  this  report  was  strictly  theoretical,  an  experimental 
analysis  should  be  made  to  verify  the  results  obtained  herein. 

2.  It  would  be  informative  to  do  both  theoretical  and  experimental 
research  on  the  effect  of  extremely  flexible  blades  upon  forward  flight. 

3.  Although  this  report  only  considered  classical  futter,  the  other 
two  types  of  flutter  -  wake  flutter  and  stall  flutter  -  should  certainly 
be  looked  at. 

4.  One  thing  that  probably  will  be  a  big  factor  in  the  construction 
of  extremely  flexible  rotor  blades  is  a  tip  mass.  Therefore,  the  effect 
of  a  tip  mass  upon  the  flutter  boundaries  of  a  rotor  blade  is  in  itself  a 
big  subject  to  consider.  Also,  a  mass  distribution,  other  than  uniform 
might  be  considered. 
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TABLE  1.   DYNAMIC  STABILITY  OF  CASE  (1) 


MODE  SHAPE 

DYNAMIC 

STABILITY 

1 

unstable 

for  '^o/fi  <  5, 

stable  elsewhere 

2 

unstable 

for  '"^/fi  <  6, 

stable  elsewhere 

3 

unstable 

for  ^°/fi  <  6, 

stable  elsewhere 

4 

unstable 

for  '^^/^  <  7, 

stable  elsewhere 
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F/GURE    2.       BLADE    TORSIONAL 
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FIGURE   23.     FLAPPIN&      MODES 
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APPENDIX  I --DETERMINANT  EXPANSION 


^0 


EXPANDED  5  THE  CHARMTERiST IC     EQUATION 
(S    OF   THE    FORM: 

Av^f  Bv^f  CV^  f  D-V'  f  Ev"  f  FVi-  6 

WHERE, 


i-!-' 


C-  -Ix^^hA^    -f-rx^gM^-Iy   JjV^^  f  ZxXa'^e  "  2Xxz  -t-J.x^•'rno^- 


t  li  M  ©  O^^   -Xlz  yy^k   T  XI2.  (^)  ^    t  tA^.^a^\     i-  Z  X2_\^  /7-1 


/ru     I     /»  I  ©  -ma  T^rt 
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E  -   -   2.Tx   Lc^^   *-  ^o  ^^  M^    /■  1a  X-z^mo  -  ^e  mc^  /\Ap   +  Lx~Lz  -K^S^I 
+  Uz.  C^)     +  ^Ae  ^^5   +  M^  'r^oi.^i.  -  -Lxi  M^  /yna.  -  l\A&rrit  -l^^t 


-rc^' 


/?^i 


■*-  Xjc  I2./V\^  ■''^e  ■/•  Xz  o^o^.^^  -  Xyi.Xz^mf  +  Xi  Me  ^(«  V'i   ^  Xli-m^  ^^ 

^llzi^fr^^<}i    ~^G^    M^     tXx^M^     f    Me^^     f  X^3    T^XC^'iVl^ 
+  X2  M4   ^?2^-Xx'^e.^'mi 

G^X/Xz-^sV-T-XAlaVi  -f  Xz^^i  -Xx^,Xi  +  XaMeS?i   -^II2^?I: 

^xx.(av 
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APPENDIX  II--ROUTH'S  ARRAY 
To  establish  stability  of  the  rotating  helicopter  blade,  with 
three  modes  of  motion,  Routh's  criteria  was  used.   In  order  to  use 
Routh's  Criteria,  Routh's  array  had  to  be  set  up.  The  form  of  the 
characteristic  equation  was: 

^^%  3VV  Cv'^^r  Dv^^t  £■^7^^f^?f  6-6  Cm) 

This  leads  to  an  array  of  the  following  form: 
.6 


v" 

A 

C 

E 

v5 

B 

D 

F 

V^ 

\ 

■^1 

G 

v^ 

^2 

h 

v^ 

N 

G 

vl 

\ 

v° 

G 

where 

\   - 

BC 

B 

AD 

h    - 

BE 

- 

AF 

^2  =  "1°  -  ^h 


\ 


B  =  A.F  -  BG 


^ 


*3  =  ^A_l^^ 
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and 


A.    =  A„B-    -   A_G 
4  ^2  2 
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APPENDIX  III^-ROTOR  BLADE  DATA 


Cornell  Aero  Lab  Blade  -^1 
NACA  0012  Airfoil 
Radius  =  48  inches 
Chord  =  3.5  inches 

Woij.;liL  "  2.i>  j)uiuiLl:i, 

2 
r  =  .000.479  slug  feet 

3 
I,  =  1.66  slug  feet 

^  =  2.46 
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APPENDIX  IV--~TERM5  OF  EQUATIONS  OF  MOTION 
r.  Flapping  and  Flap-Bending  Equations 

1 .  3  and  g  terms . 

Tkese  terms  arise  from  the  z  -  component  of  the  centrifugal  force 
at  each  blade  element.  A  virtual  deflection  of  63  or  6g  results  in 
energy  being  expended  which  results  in  a  force  for  the  whole  blade  of 

.a^Tj^^  CIVI) 

2.  e  term. 

Due  to  a  twist  of  the  blade,  the  centrifugal  force  of  a  blade  element, 

2 

Mrfi  dr,  will  act  with  a  moment  arm,  -  it,  X  dr9,  resulting  in  a  force  of 

-j'^Tr.a^q^X^d^  ©  CIV3) 

3.  B  and  g  term. 

At  each  blade  element  there  is  an  intertial  force  of 

which  results  in  a  force  of 

-hi 

M2.3  CIVS) 

4.  0  term. 

Due  to  a  pitching  acceleration  there  is  an  inertial  effect  which  re- 
sults in  the  following  force. 

-J^  nrnn.  f]^^XJ.  dy    e  CIV61 
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5.   0  term,  aerodynamic. 
This  term  is  due  to  a  twist  angle,  9,  which  causes  energy  to  be 
transmitted  by  the  airstream  to  the  blade  giving  a  force  of 


CrV7) 
and 


It  ^''^o  S 


6.   g  and  g  terms. 
The  flapping  velocity  at  each,  blade  element  causes  a  relative  change 
in  the  angle  of  attack  which  results  in  the  following  forces: 

Xi^.  '^  Crv9} 

-r2.^cj  -J  Civio) 

^^^^^  Civil) 

^^2,^  Crvi2) 

ir.   Pitching  Equation 

1.  6  term. 

Due  to  a  virtual  twist  of  the  blade,  the  z  -  component  of  the  cen- 
trifugal force  of  a  blade  element  causes  a  force  of 

'jf-rr.a^'^^K^ir  ^  CIV13) 

2.  9  terms. 

This  term  arises  from  the  "tennis  racket"  effect.   Each  partical  of 
mass,  dm,  of  a  blade  section  is  acted  upon  -  due  to  a  twist  deflection  -  by 
a  component  of  the  centrifugal  force  in  the  plane  of  the  section.  The  re- 
sulting force  is 

[r^n^  t  I^u^'.^Jg  CIV14) 
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3,  g  term. 

Due  to  a  virtual  deflection  of  a  blade  section,  the  inertia  force, 
mrgdr,  results  in  a  force  of 

4,  6  term. 

Tliis   force  arises  from  the  inertia  force  du  e  to  a  twist  of  the  blade. 

X  ii^'e  CIV16) 

5.  I  term. 

Tliis  aerodynamic  term  is  the  result  of  a  couple  produced  by  the 
pitching  acceleration  of  the  rotor  blade. 

M-  I,  a  ^  CIV17) 

6.  0  term. 

This  aerodynamic  term  is  the  result  of  the  energy  transmitted  to  the 
blade  by  the  lift  forces  moving  through  a  virtual  displacement. 

Me -^4-^^  CIV18) 
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